Genetic modifiers of rare variants in monogenic developmental disorder loci.

Rare damaging variants in a large number of genes are known to cause monogenic developmental disorders (DDs) and have also been shown to cause milder subclinical phenotypes in population cohorts. Here, we show that carrying multiple (2-5) rare damaging variants across 599 dominant DD genes has an additive adverse effect on numerous cognitive and socioeconomic traits in UK Biobank, which can be partially counterbalanced by a higher educational attainment polygenic score (EA-PGS). Phenotypic deviators from expected EA-PGS could be partly explained by the enrichment or depletion of rare DD variants. Among carriers of rare DD variants, those with a DD-related clinical diagnosis had a substantially lower EA-PGS and more severe phenotype than those without a clinical diagnosis. Our results suggest that the overall burden of both rare and common variants can modify the expressivity of a phenotype, which may then influence whether an individual reaches the threshold for clinical disease.

Rare variant associations with birth weight identify genes involved in adipose tissue regulation, placental function and insulin-like growth factor signalling.

Investigating the genetic factors influencing human birth weight may lead to biological insights into fetal growth and long-term health. Genome-wide association studies of birth weight have highlighted associated variants in more than 200 regions of the genome, but the causal genes are mostly unknown. Rare genetic variants with robust evidence of association are more likely to point to causal genes, but to date, only a few rare variants are known to influence birth weight. We aimed to identify genes that harbour rare variants that impact birth weight when carried by either the fetus or the mother, by analysing whole exome sequence data in UK Biobank participants. We annotated rare (minor allele frequency <0.1%) protein-truncating or high impact missense variants on whole exome sequence data in up to 234,675 participants with data on their own birth weight (fetal variants), and up to 181,883 mothers who reported the birth weight of their first child (maternal variants). Variants within each gene were collapsed to perform gene burden tests and for each associated gene, we compared the observed fetal and maternal effects. We identified 8 genes with evidence of rare fetal variant effects on birth weight, of which 2 also showed maternal effects. One additional gene showed evidence of maternal effects only. We observed 10/11 directionally concordant associations in an independent sample of up to 45,622 individuals (sign test P=0.01). Of the genes identified, IGF1R and PAPPA2 (fetal and maternal-acting) have known roles in insulin-like growth factor bioavailability and signalling. PPARG, INHBE and ACVR1C (all fetal-acting) have known roles in adipose tissue regulation and rare variants in the latter two also showed associations with favourable adiposity patterns in adults. We highlight the dual role of PPARG in both adipocyte differentiation and placental angiogenesis. NOS3, NRK, and ADAMTS8 (fetal and maternal-acting) have been implicated in both placental function and hypertension. Analysis of rare coding variants has identified regulators of fetal adipose tissue and fetoplacental angiogenesis as determinants of birth weight, as well as further evidence for the role of insulin-like growth factors.

Hyperglycaemia is a causal risk factor for upper limb pathologies.

BACKGROUND: Diabetes (regardless of type) and obesity are associated with a range of musculoskeletal disorders. The causal mechanisms driving these associations are unknown for many upper limb pathologies. We used genetic techniques to test the causal link between glycemia, obesity and musculoskeletal conditions. METHODS: In the UK Biobank's unrelated European cohort (N = 379 708) we performed mendelian randomisation (MR) analyses to test for a causal effect of long-term high glycaemia and adiposity on four musculoskeletal pathologies: frozen shoulder, Dupuytren's disease, carpal tunnel syndrome and trigger finger. We also performed single-gene MR using rare variants in the GCK gene. RESULTS: Using MR, we found evidence that long-term high glycaemia has a causal role in the aetiology of upper limb conditions. A 10-mmol/mol increase in genetically predicted haemoglobin A1C (HbA1c) was associated with frozen shoulder: odds ratio (OR) = 1.50 [95% confidence interval (CI), 1.20-1.88], Dupuytren's disease: OR = 1.17 (95% CI, 1.01-1.35), trigger finger: OR = 1.30 (95% CI, 1.09-1.55) and carpal tunnel syndrome: OR = 1.20 (95% CI, 1.09-1.33). Carriers of GCK mutations have increased odds of frozen shoulder: OR = 7.16 (95% CI, 2.93-17.51) and carpal tunnel syndrome: OR = 2.86 (95% CI, 1.50-5.44) but not Dupuytren's disease or trigger finger. We found evidence that an increase in genetically predicted body mass index (BMI) of 5 kg/m2 was associated with carpal tunnel syndrome: OR = 1.13 (95% CI, 1.10-1.16) and associated negatively with Dupuytren's disease: OR = 0.94 (95% CI, 0.90-0.98), but no evidence of association with frozen shoulder or trigger finger. Trigger finger (OR 1.96 (95% CI, 1.42-2.69) P = 3.6e-05) and carpal tunnel syndrome [OR 1.63 (95% CI, 1.36-1.95) P = 8.5e-08] are associated with genetically predicted unfavourable adiposity increase of one standard deviation of body fat. CONCLUSIONS: Our study consistently demonstrates a causal role of long-term high glycaemia in the aetiology of upper limb musculoskeletal conditions. Clinicians treating diabetes patients should be aware of these complications in clinic, specifically those managing the care of GCK mutation carriers. Upper limb musculoskeletal conditions should be considered diabetes complications.

Iron and risk of dementia: Mendelian randomisation analysis in UK Biobank.

BACKGROUND: Brain iron deposition is common in dementia, but whether serum iron is a causal risk factor is unknown. We aimed to determine whether genetic predisposition to higher serum iron status biomarkers increased risk of dementia and atrophy of grey matter. METHODS: We analysed UK Biobank participants clustered into European (N=451284), African (N=7477) and South Asian (N=9570) groups by genetic similarity to the 1000 genomes project. Using Mendelian randomisation methods, we estimated the association between genetically predicted serum iron (transferrin saturation [TSAT] and ferritin), grey matter volume and genetic liability to clinically defined dementia (including Alzheimer's disease [AD], non-AD dementia, and vascular dementia) from hospital and primary care records. We also performed time-to-event (competing risks) analysis of the TSAT polygenic score on risk of clinically defined non-AD dementia. RESULTS: In Europeans, higher genetically predicted TSAT increased genetic liability to dementia (Odds Ratio [OR]: 1.15, 95% Confidence Intervals [CI] 1.04 to 1.26, p=0.0051), non-AD dementia (OR: 1.27, 95% CI 1.12 to 1.45, p=0.00018) and vascular dementia (OR: 1.37, 95% CI 1.12 to 1.69, p=0.0023), but not AD (OR: 1.00, 95% CI 0.86 to 1.15, p=0.97). Higher TSAT was also associated with increased risk of non-AD dementia in participants of African, but not South Asian groups. In survival analysis using a TSAT polygenic score, the effect was independent of apolipoprotein-E ε4 genotype (with adjustment subdistribution Hazard Ratio: 1.74, 95% CI 1.33 to 2.28, p=0.00006). Genetically predicted TSAT was associated with lower grey matter volume in caudate, putamen and thalamus, and not in other areas of interest. DISCUSSION: Genetic evidence supports a causal relationship between higher TSAT and risk of clinically defined non-AD and vascular dementia, in European and African groups. This association appears to be independent of apolipoprotein-E ε4.

Does physical activity moderate the association between shorter leukocyte telomere length and incident coronary heart disease? Data from 54,180 UK Biobank participants.

Telomere shortening is a biological aging hallmark. The effect of short telomere length may be targeted by increased physical activity to reduce the risk of multiple aging-related diseases, including coronary heart disease (CHD). The objective was to assess the moderation effect of accelerometer-based physical activity (aPA) on the association between shorter leukocyte telomere length (LTL) relatively in the population sample and incident CHD. Data were from the UK Biobank participants with well-calibrated accelerometer data for at least 6.5 days (n = 54,180). Relative mean LTL at baseline (5-6 years prior to aPA assessment) was measured in T/S ratio, using a multiplex quantitative polymerase chain reaction (qPCR) technology, by comparing the amount of the telomere amplification product (T) to that of a single-copy gene (S). aPA measures included total number of events (at least 10-s continued physical activity > 32 milligravities [mg]), total volume, mean duration, mean intensity, and peak intensity of all events. LTL, aPA measures, and their interactions were associated with incident CHD (mean follow-up 6.8 years) using Cox proportional hazards models adjusting for covariates. Longer LTL (relative to the sample distribution) was associated with reduced incidence of CHD (adjusted hazard ratio [aHR] = 0.94 per standard deviation [SD] increase in LTL, [95% CI, 0.90 to 0.99], P = .010). Incidence of CHD was reduced by higher total volume of aPA (aHR = 0.82 per SD increase in LTL, [95% CI, 0.71 to 0.95], P = .010) but increased by higher total number of events (aHR = 1.11 per SD increase in LTL, [95% CI, 1.02 to 1.21], P = .020) after controlling for other aPA measures and covariates. However, none of the interactions between LTL and aPA measures was statistically significant (P = .171).

Effects of physical activity and sedentary time on depression, anxiety and well-being: a bidirectional Mendelian randomisation study.

BACKGROUND: Mental health conditions represent one of the major groups of non-transmissible diseases. Physical activity (PA) and sedentary time (ST) have been shown to affect mental health outcomes in opposite directions. In this study, we use accelerometery-derived measures of PA and ST from the UK Biobank (UKB) and depression, anxiety and well-being data from the UKB mental health questionnaire as well as published summary statistics to explore the causal associations between these phenotypes. METHODS: We used MRlap to test if objectively measured PA and ST associate with mental health outcomes using UKB data and summary statistics from published genome-wide association studies. We also tested for bidirectional associations. We performed sex stratified as well as sensitivity analyses. RESULTS: Genetically instrumented higher PA was associated with lower odds of depression (OR = 0.92; 95% CI: 0.88, 0.97) and depression severity (beta = - 0.11; 95% CI: - 0.18, - 0.04), Genetically instrumented higher ST was associated higher odds of anxiety (OR = 2.59; 95% CI: 1.10, 4.60). PA was associated with higher well-being (beta = 0.11, 95% CI: 0.04; 0.18) and ST with lower well-being (beta = - 0.18; 95% CI: - 0.32, - 0.03). Similar findings were observed when stratifying by sex. There was evidence for a bidirectional relationship, with higher genetic liability to depression associated with lower PA (beta = - 0.25, 95% CI: - 0.42; - 0.08) and higher well-being associated with higher PA (beta = 0.15; 95% CI: 0.05, 0.25). CONCLUSIONS: We have demonstrated the bidirectional effects of both PA and ST on a range of mental health outcomes using objectively measured predictors and MR methods for causal inference. Our findings support a causal role for PA and ST in the development of mental health problems and in affecting well-being.

Insights into the genetics of menopausal vasomotor symptoms: genome-wide analyses of routinely-collected primary care health records.

BACKGROUND: Vasomotor symptoms (VMS) can often significantly impact women's quality of life at menopause. In vivo studies have shown that increased neurokinin B (NKB) / neurokinin 3 receptor (NK3R) signalling contributes to VMS, with previous genetic studies implicating the TACR3 gene locus that encodes NK3R. Large-scale genomic analyses offer the possibility of biological insights but few such studies have collected data on VMS, while proxy phenotypes such as hormone replacement therapy (HRT) use are likely to be affected by changes in clinical practice. We investigated the genetic basis of VMS by analysing routinely-collected health records. METHODS: We performed a GWAS of VMS derived from linked primary-care records and cross-sectional self-reported HRT use in up to 153,152 women from UK Biobank, a population-based cohort. In a subset of this cohort (n = 39,356), we analysed exome-sequencing data to test the association with VMS of rare deleterious genetic variants. Finally, we used Mendelian randomisation analysis to investigate the reasons for HRT use over time. RESULTS: Our GWAS of health-records derived VMS identified a genetic signal near TACR3 associated with a lower risk of VMS (OR=0.76 (95% CI 0.72,0.80) per A allele, P=3.7x10-27), which was consistent with previous studies, validating this approach. Conditional analyses demonstrated independence of genetic signals for puberty timing and VMS at the TACR3 locus, including a rare variant predicted to reduce functional NK3R levels that was associated with later menarche (P = 5 × 10-9) but showed no association with VMS (P = 0.6). Younger menopause age was causally-associated with greater HRT use before 2002 but not after. CONCLUSIONS: We provide support for TACR3 in the genetic basis of VMS but unexpectedly find that rare genomic variants predicted to lower NK3R levels did not modify VMS, despite the proven efficacy of NK3R antagonists. Using genomics we demonstrate changes in genetic associations with HRT use over time, arising from a change in clinical practice since the early 2000s, which is likely to reflect a switch from preventing post-menopausal complications in women with earlier menopause to primarily treating VMS. Our study demonstrates that integrating routinely-collected primary care health records and genomic data offers great potential for exploring the genetic basis of symptoms.

Identification and analysis of individuals who deviate from their genetically-predicted phenotype.

Findings from genome-wide association studies have facilitated the generation of genetic predictors for many common human phenotypes. Stratifying individuals misaligned to a genetic predictor based on common variants may be important for follow-up studies that aim to identify alternative causal factors. Using genome-wide imputed genetic data, we aimed to classify 158,951 unrelated individuals from the UK Biobank as either concordant or deviating from two well-measured phenotypes. We first applied our methods to standing height: our primary analysis classified 244 individuals (0.15%) as misaligned to their genetically predicted height. We show that these individuals are enriched for self-reporting being shorter or taller than average at age 10, diagnosed congenital malformations, and rare loss-of-function variants in genes previously catalogued as causal for growth disorders. Secondly, we apply our methods to LDL cholesterol (LDL-C). We classified 156 (0.12%) individuals as misaligned to their genetically predicted LDL-C and show that these individuals were enriched for both clinically actionable cardiovascular risk factors and rare genetic variants in genes previously shown to be involved in metabolic processes. Individuals whose LDL-C was higher than expected based on the genetic predictor were also at higher risk of developing coronary artery disease and type-two diabetes, even after adjustment for measured LDL-C, BMI and age, suggesting upward deviation from genetically predicted LDL-C is indicative of generally poor health. Our results remained broadly consistent when performing sensitivity analysis based on a variety of parametric and non-parametric methods to define individuals deviating from polygenic expectation. Our analyses demonstrate the potential importance of quantitatively identifying individuals for further follow-up based on deviation from genetic predictions.

Penetrance of pathogenic genetic variants associated with premature ovarian insufficiency.

Premature ovarian insufficiency (POI) affects 1% of women and is a leading cause of infertility. It is often considered to be a monogenic disorder, with pathogenic variants in ~100 genes described in the literature. We sought to systematically evaluate the penetrance of variants in these genes using exome sequence data in 104,733 women from the UK Biobank, 2,231 (1.14%) of whom reported at natural menopause under the age of 40 years. We found limited evidence to support any previously reported autosomal dominant effect. For nearly all heterozygous effects on previously reported POI genes, we ruled out even modest penetrance, with 99.9% (13,699 out of 13,708) of all protein-truncating variants found in reproductively healthy women. We found evidence of haploinsufficiency effects in several genes, including TWNK (1.54 years earlier menopause, P = 1.59 × 10-6) and SOHLH2 (3.48 years earlier menopause, P = 1.03 × 10-4). Collectively, our results suggest that, for the vast majority of women, POI is not caused by autosomal dominant variants either in genes previously reported or currently evaluated in clinical diagnostic panels. Our findings, plus previous studies, suggest that most POI cases are likely oligogenic or polygenic in nature, which has important implications for future clinical genetic studies, and genetic counseling for families affected by POI.

Identification and analysis of individuals who deviate from their genetically-predicted phenotype.

Findings from genome-wide association studies have facilitated the generation of genetic predictors for many common human phenotypes. Stratifying individuals misaligned to a genetic predictor based on common variants may be important for follow-up studies that aim to identify alternative causal factors. Using genome-wide imputed genetic data, we aimed to classify 158,951 unrelated individuals from the UK Biobank as either concordant or deviating from two well-measured phenotypes. We first applied our methods to standing height: our primary analysis classified 244 individuals (0.15%) as misaligned to their genetically predicted height. We show that these individuals are enriched for self-reporting being shorter or taller than average at age 10, diagnosed congenital malformations, and rare loss-of-function variants in genes previously catalogued as causal for growth disorders. Secondly, we apply our methods to LDL cholesterol. We classified 156 (0.12%) individuals as misaligned to their genetically predicted LDL cholesterol and show that these individuals were enriched for both clinically actionable cardiovascular risk factors and rare genetic variants in genes previously shown to be involved in metabolic processes. Individuals whose LDL-C was higher than expected based on the genetic predictor were also at higher risk of developing coronary artery disease and type-two diabetes, even after adjustment for measured LDL-C, BMI and age, suggesting upward deviation from genetically predicted LDL-C is indicative of generally poor health. Our results remained broadly consistent when performing sensitivity analysis based on a variety of parametric and non-parametric methods to define individuals deviating from polygenic expectation. Our analyses demonstrate the potential importance of quantitatively identifying individuals for further follow-up based on deviation from genetic predictions. Author Summary: Human genetics is becoming increasingly useful to help predict human traits across a population owing to findings from large-scale genetic association studies and advances in the power of genetic predictors. This provides an opportunity to potentially identify individuals that deviate from genetic predictions for a common phenotype under investigation. For example, an individual may be genetically predicted to be tall, but be shorter than expected. It is potentially important to identify individuals who deviate from genetic predictions as this can facilitate further follow-up to assess likely causes. Using 158,951 unrelated individuals from the UK Biobank, with height and LDL cholesterol, as exemplar traits, we demonstrate that approximately 0.15% & 0.12% of individuals deviate from their genetically predicted phenotypes respectively. We observed these individuals to be enriched for a range of rare clinical diagnoses, as well as rare genetic factors that may be causal. Our analyses also demonstrate several methods for detecting individuals who deviate from genetic predictions that can be applied to a range of continuous human phenotypes.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk.

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

Recurrent 17q12 microduplications contribute to renal disease but not diabetes.

BACKGROUND: 17q12 microdeletion and microduplication syndromes present as overlapping, multisystem disorders. We assessed the disease phenotypes of individuals with 17q12 CNV in a population-based cohort. METHODS: We investigated 17q12 CNV using microarray data from 450 993 individuals in the UK Biobank and calculated disease status associations for diabetes, liver and renal function, neurological and psychiatric traits. RESULTS: We identified 11 17q12 microdeletions and 106 microduplications. Microdeletions were strongly associated with diabetes (p=2×10-7) but microduplications were not. Estimated glomerular filtration rate (eGFR mL/min/1.73 m2) was consistently lower in individuals with microdeletions (p=3×10-12) and microduplications (p=6×10-25). Similarly, eGFR <60, including end-stage renal disease, was associated with microdeletions (p=2×10-9, p<0.003) and microduplications (p=1×10-9, p=0.009), respectively, highlighting sometimes substantially reduced renal function in each. Microduplications were associated with decreased fluid intelligence (p=3×10-4). SNP association analysis in the 17q12 region implicated changes to HNF1B as causing decreased eGFR (NC_000017.11:g.37741642T>G, rs12601991, p=4×10-21) and diabetes (NC_000017.11:g.37741165C>T, rs7501939, p=6×10-17). A second locus within the region was also associated with fluid intelligence (NC_000017.11:g.36593168T>C, rs1005552, p=6×10-9) and decreased eGFR (NC_000017.11:g.36558947T>C, rs12150665, p=4×10-15). CONCLUSION: We demonstrate 17q12 microdeletions but not microduplications are associated with diabetes in a population-based cohort, likely caused by HNF1B haploinsufficiency. We show that both 17q12 microdeletions and microduplications are associated with renal disease, and multiple genes within the region likely contribute to renal and neurocognitive phenotypes.

Genetically proxied therapeutic prolyl-hydroxylase inhibition and cardiovascular risk.

Prolyl hydroxylase (PHD) inhibitors are in clinical development for anaemia in chronic kidney disease. Epidemiological studies have reported conflicting results regarding safety of long-term therapeutic haemoglobin (Hgb) rises through PHD inhibition on risk of cardiovascular disease. Genetic variation in genes encoding PHDs can be used as partial proxies to investigate the potential effects of long-term Hgb rises. We used Mendelian randomization to investigate the effect of long-term Hgb level rises through genetically proxied PHD inhibition on coronary artery disease (CAD: 60 801 cases; 123 504 controls), myocardial infarction (MI: 42 561 cases; 123 504 controls) or stroke (40 585 cases; 406 111 controls). To further characterize long-term effects of Hgb level rises, we performed a phenome-wide association study (PheWAS) in up to 451 099 UK Biobank individuals. Genetically proxied therapeutic PHD inhibition, equivalent to a 1.00 g/dl increase in Hgb levels, was not associated (at P < 0.05) with increased odds of CAD; odd ratio (OR) [95% confidence intervals (CI)] = 1.06 (0.84, 1.35), MI [OR (95% CI) = 1.02 (0.79, 1.33)] or stroke [OR (95% CI) = 0.91 (0.66, 1.24)]. PheWAS revealed associations with blood related phenotypes consistent with EGLN's role, relevant kidney- and liver-related biomarkers like estimated glomerular filtration rate and microalbuminuria, and non-alcoholic fatty liver disease (Bonferroni-adjusted P < 5.42E-05) but these were not clinically meaningful. These findings suggest that long-term alterations in Hgb through PHD inhibition are unlikely to substantially increase cardiovascular disease risk; using large disease genome-wide association study data, we could exclude ORs of 1.35 for cardiovascular risk with a 1.00 g/dl increase in Hgb.

Prevalence of Fabry disease-causing variants in the UK Biobank.

BACKGROUND: Fabry disease is an X-linked lysosomal storage disorder resulting from deficiency of the alpha-galactosidase A enzyme leading to accumulation of globotriaosylceramide in multiple organ sites with prominent cardiovascular and renal involvement. Global prevalence estimates of Fabry disease based on clinical ascertainment range from 1 in 40 000 to 1 in 170 000. We aimed to determine the prevalence of Fabry disease-causing variants in UK Biobank. METHODS: We sought GLA gene variants in exome sequencing data from 200 643 individuals from UK Biobank. We used ACMG/AMP guidelines (American College of Medical Genetics/Association for Molecular Pathology) to classify pathogenicity and compared baseline biomarker data, hospital ICD-10 (International Classification of Diseases version-10) codes, general practitioner records and self-reported health data with those without pathogenic variants. RESULTS: We identified 81 GLA coding variants. We identified eight likely pathogenic variants on the basis of being rare (<1/10 000 individuals) and either previously reported to cause Fabry disease, or being protein-truncating variants. Thirty-six individuals carried one of these variants. In the UK Biobank, the prevalence of likely pathogenic Fabry disease-causing variants is 1/5732 for late-onset disease-causing variants and 1/200 643 for variants causing classic Fabry disease. CONCLUSION: Fabry disease-causing GLA variants are more prevalent in an unselected population sample than the reported prevalence of Fabry disease. These are overwhelmingly variants associated with later onset. It is possible the prevalence of later-onset Fabry disease exceeds current estimates.

Reduced penetrance of MODY-associated HNF1A/HNF4A variants but not GCK variants in clinically unselected cohorts.

The true prevalence and penetrance of monogenic disease variants are often not known because of clinical-referral ascertainment bias. We comprehensively assess the penetrance and prevalence of pathogenic variants in HNF1A, HNF4A, and GCK that account for >80% of monogenic diabetes. We analyzed clinical and genetic data from 1,742 clinically referred probands, 2,194 family members, clinically unselected individuals from a US health system-based cohort (n = 132,194), and a UK population-based cohort (n = 198,748). We show that one in 1,500 individuals harbor a pathogenic variant in one of these genes. The penetrance of diabetes for HNF1A and HNF4A pathogenic variants was substantially lower in the clinically unselected individuals compared to clinically referred probands and was dependent on the setting (32% in the population, 49% in the health system cohort, 86% in a family member, and 98% in probands for HNF1A). The relative risk of diabetes was similar across the clinically unselected cohorts highlighting the role of environment/other genetic factors. Surprisingly, the penetrance of pathogenic GCK variants was similar across all cohorts (89%-97%). We highlight that pathogenic variants in HNF1A, HNF4A, and GCK are not ultra-rare in the population. For HNF1A and HNF4A, we need to tailor genetic interpretation and counseling based on the setting in which a pathogenic monogenic variant was identified. GCK is an exception with near-complete penetrance in all settings. This along with the clinical implication of diagnosis makes it an excellent candidate for the American College of Medical Genetics secondary gene list.

The impact of Mendelian sleep and circadian genetic variants in a population setting.

Rare variants in ten genes have been reported to cause Mendelian sleep conditions characterised by extreme sleep duration or timing. These include familial natural short sleep (ADRB1, DEC2/BHLHE41, GRM1 and NPSR1), advanced sleep phase (PER2, PER3, CRY2, CSNK1D and TIMELESS) and delayed sleep phase (CRY1). The association of variants in these genes with extreme sleep conditions were usually based on clinically ascertained families, and their effects when identified in the population are unknown. We aimed to determine the effects of these variants on sleep traits in large population-based cohorts. We performed genetic association analysis of variants previously reported to be causal for Mendelian sleep and circadian conditions. Analyses were performed using 191,929 individuals with data on sleep and whole-exome or genome-sequence data from 4 population-based studies: UK Biobank, FINRISK, Health-2000-2001, and the Multi-Ethnic Study of Atherosclerosis (MESA). We identified sleep disorders from self-report, hospital and primary care data. We estimated sleep duration and timing measures from self-report and accelerometery data. We identified carriers for 10 out of 12 previously reported pathogenic variants for 8 of the 10 genes. They ranged in frequency from 1 individual with the variant in CSNK1D to 1,574 individuals with a reported variant in the PER3 gene in the UK Biobank. No carriers for variants reported in NPSR1 or PER2 were identified. We found no association between variants analyzed and extreme sleep or circadian phenotypes. Using sleep timing as a proxy measure for sleep phase, only PER3 and CRY1 variants demonstrated association with earlier and later sleep timing, respectively; however, the magnitude of effect was smaller than previously reported (sleep midpoint ~7 mins earlier and ~5 mins later, respectively). We also performed burden tests of protein truncating (PTVs) or rare missense variants for the 10 genes. Only PTVs in PER2 and PER3 were associated with a relevant trait (for example, 64 individuals with a PTV in PER2 had an odds ratio of 4.4 for being "definitely a morning person", P = 4x10-8; and had a 57-minute earlier midpoint sleep, P = 5x10-7). Our results indicate that previously reported variants for Mendelian sleep and circadian conditions are often not highly penetrant when ascertained incidentally from the general population.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies.

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Detection and characterization of male sex chromosome abnormalities in the UK Biobank study.

PURPOSE: The study aimed to systematically ascertain male sex chromosome abnormalities, 47,XXY (Klinefelter syndrome [KS]) and 47,XYY, and characterize their risks of adverse health outcomes. METHODS: We analyzed genotyping array or exome sequence data in 207,067 men of European ancestry aged 40 to 70 years from the UK Biobank and related these to extensive routine health record data. RESULTS: Only 49 of 213 (23%) of men whom we identified with KS and only 1 of 143 (0.7%) with 47,XYY had a diagnosis of abnormal karyotype on their medical records or self-report. We observed expected associations for KS with reproductive dysfunction (late puberty: risk ratio [RR] = 2.7; childlessness: RR = 4.2; testosterone concentration: RR = -3.8 nmol/L, all P < 2 × 10-8), whereas XYY men appeared to have normal reproductive function. Despite this difference, we identified several higher disease risks shared across both KS and 47,XYY, including type 2 diabetes (RR = 3.0 and 2.6, respectively), venous thrombosis (RR = 6.4 and 7.4, respectively), pulmonary embolism (RR = 3.3 and 3.7, respectively), and chronic obstructive pulmonary disease (RR = 4.4 and 4.6, respectively) (all P < 7 × 10-6). CONCLUSION: KS and 47,XYY were mostly unrecognized but conferred substantially higher risks for metabolic, vascular, and respiratory diseases, which were only partially explained by higher levels of body mass index, deprivation, and smoking.

Simulated distributions from negative experiments highlight the importance of the body mass index distribution in explaining depression-body mass index genetic risk score interactions.

BACKGROUND: Depression and obesity are complex global health problems. Recent studies suggest that a genetic predisposition to obesity might be accentuated in people with depression, but these analyses are prone to bias. Here, we tested the hypothesis that depression accentuates genetic susceptibility to obesity and applied negative control experiments to test whether any observed interactions were real or driven by confounding and statistical biases. METHODS: We used data from up to 378 000 Europeans in UK Biobank, a 73 variant body mass index (BMI) genetic risk score, two depression measures [depression symptoms (DS), major depression (MD)] and an antidepressant usage variable available. We tested whether (i) depression and (ii) antidepressant treatment accentuated genetic susceptibility to obesity. Finally, we performed negative control experiments by sampling individuals at random so that they had BMI distributions identical to depression cases and controls. RESULTS: Depression was associated with an accentuation of an individual's genetic risk of obesity with evidence of interactions for both DS and MD (Pinteraction = 7 × 10-4 and 7 × 10-5 respectively). Antidepressant usage within DS cases accentuated genetic obesity risk (Pinteraction = 9 × 10-4), but not for MD (Pinteraction = 0.13). Negative control experiments suggested that the observed interactions for MD (empirical-P = 0.067) may be driven by statistical biases or confounding factors but were not possible with the larger DS groups. Antidepressant usage interaction also appears to be driven by statistical artefacts (empirical-P = 0.510 using MD and 0.162 using DS). CONCLUSION: We have highlighted the importance of running negative experiments to confirm putative interactions in gene-environment studies. We provide some tentative evidence that depression accentuates an individual's genetic susceptibility to higher BMI but demonstrated that the BMI distributions within cases and controls might drive these interactions.